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Abstract 

The Holstein-Primakoff transformation is generalized to develop a quan- 
tum flavor wave theory for spin systems with orbital degeneracy. Elemen- 
tary excitations of ordered ground states consist of spin, orbital , and spin- 
orbital waves. Spin and spin-orbital waves couple to each other due to orbital 
anisotropy and Hund's rule, resulting in new modes observable by inelastic 
neutron scattering. In the SU{4) limit, flavor waves are dispersionless along 
one or more directions, and give rise to quantum fluctuations of reduced di- 
mensionality. 
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The Hubbard model with a Hilbert space of only one atomic orbital per site has been a 
popular model used to study strongly correlated electronic systems. When the correlation U 
is sufficiently large, the system undergoes a transition into a Mott insulator. In the large U 
limit, the model reduces to the familiar S — 1/2 Heisenberg antiferromagnetic Hamiltonian 
(HAH). In spite of the success and popularity of the single orbital Hubbard model, it cannot, 
however, explain the magnetic behavior of many Mott insulators, including ferromagnetism, 
magnetic ordering patterns, such as that in V2O3 [PHI, which would require unrealistically 
large long-ranged coupling if explained by HAH, and paramagnetic behavior on lattices 
where the HAH is known to exhibit long-range order (LRO) [|J. On the other hand, such 
behavior can be easily understood if one allows for orbital degeneracy and enlarges the 
Hilbert space on each site 0-01 . Recently, direct evidence for orbital degrees of freedom 
and ordering in some magnetic systems have been obtained by form factor analysis of x-ray 
diffraction || . 

In the case of spin systems with two-fold orbital degeneracy, the Hamiltonian is quadratic 
in operators S™ , the generators of SU(4) Lie algebra. Indeed the Hamiltonian may have 
SU(4) symmetry [|T|,|TU|, rather than just the spin SU(2) symmetry. In many 2D lattices 
such as the square or triangular lattice with nearest neighbor (n.n) coupling, the 5?7(4) 
antiferromagnetic (AF) Hamiltonian has no LRO even at T = 0, as indicated by mean field 
theory, variational calculations J7|, and recent Monte Carlo simulations 0. LRO may be 
stabilized in higher dimensions. In real systems, the 577(4) symmetry is at best approximate 
due to orbital anisotropy and Hund's rule. With sufficient deviation from the 577(4) limit, 
LRO may be attained even in 2D. 

In this paper we investigate the elementary excitations of the orbital-spin systems as- 
suming the ground state has broken symmetry. In effect, we are studying for orbital-spin 
systems the equivalence of the familiar spin waves, henceforth called flavor waves. For spin 
only systems, the Holstein-Primakoff transformation (HPT) maps spin operators into boson 
operators, and the linearized spin wave theory is equivalent to a non-interacting boson prob- 
lem. For orbital-spin systems, we show that the 577(4) algebra can be exactly reproduced 
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by a generalized HPT involving three bosons (N — 1 bosons for SU(N)). A quantum flavor 
wave theory is then developed. We find that in the 577(4) limit, even though the underlying 
lattice and ordering pattern is two-dimensional (2D), the flavor wave excitations can be 
one-dimensional (dispersionless along one direction) or localized (dispersionless in all direc- 
tions). Quantum fluctuations of these excitations give rise to disordering effect of reduced 
dimensionality which provides further support for the lack of LRO f7||| . The excitations in 
the S77(4) limit can be characterized as pure spin waves (spin rotation only), pure orbital 
waves (orbital rotation only), and pure spin-orbital waves (simultaneous spin and orbital 
rotation). The anisotropy and the Hund's rule break the SU(4) symmetry, and the spin 
and spin-orbital waves are in general mixed. We use a simplified model relevant to V2O3 to 
illustrate the mixing, and predict new modes observable in neutron scattering experiments. 

With two orbital degrees of freedom, we can define orbital operator f which acts on 
orbital states in the same way spin operator s on spin states. The Hilbert space on each site 
consists of 4 basis states, which we choose as | s z ,r z ). We label them as 

11 1 1 

I 1) -I 2'2 ) ' 1 2 > -| ~2' 2'' 

|3>=|i,«>,|4>=|4-i>. (1) 

These basis states form a fundamental representation of SU(4). The conventional S77(4) 
generators acts on the basis state | I) according to \ I) = 8 nj i \ m), and satisfies 
E m S™ = 1 and (S^y = S™. The Lie algebra is given by 

s a ,T/3, and the product s a Tp can all be expressed as linear combinations of the Sl^ J7|. For 
instance, 2s z = Si — Sf + S| — = Sf + S^, and similarly for the orbital operators. 

Ignoring anisotropy and Hund's rule, the Hubbard model with double orbital degeneracy in 
the large U limit with 1 electron per site (1/4 filling) gives rise to the effective Hamiltonian 

H = J2^S^(i)S^(j), (3) 
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with repeated indices n, m summed. H is clearly invariant under global SU(4) transfor- 
mation. We have argued J7J that the ground state of H on many 2D lattices should be a 
flavor liquid. It is however useful to study the possible excitations assuming the ground 
state has spontaneous symmetry breaking. This is because i) while the ground state for the 
fundamental representation is disordered, it may have LRO for higher representations ii) 
such a study can provide information about the stability of certain classical ground states 
against quantum fluctuations; and hi) LRO can exist in 3D and/or away from the 577(4) 
limit; and this analysis will allow insight as to what aspects of those excitations are due to 
proximity to ££7(4) symmetry. 

Just as in the £77(2) case, it is useful to generalize the s = 1/2 problem to general s, 
we can also consider representations other than the fundamental representation for ££7(4). 
In particular, we generalize it to representations denoted by Young tableaux with a single 
row but arbitrary columns M. In the limit M — > oo , the non-commutativity between the 
£^ can be ignored, and the operators become c-numbers. For general M, the Lie algebra 
(Q) can be exactly reproduced by a generalized HPT using 3-bosons b™ at each site, with 
n 7^ m. The vacuum of these bosons is the state with S™ = M, where the choice of m is 
formally arbitrary, but in practice is taken to be the ordering "direction" of the classical 
ground state. The generalized HPT is defined as (for n,l ^ m): 

s™ = b™\lM-j2bfbr, 

S l n = b™X- (4) 

Eq. (0) enables us to carry out the linear flavor wave theory by expanding the Hamiltonian 
in power of 1/M to quadratic order in b and/or and then set M — 1 for the present 
physical system. 

First consider the "ferromagnetic" (FM) case, < in @. The classical ground state, 
which is also the exact ground state, is given by S%(i) = M5 n l . Expansion of (§) leads to 
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the flavor wave Hamiltonian (omitting the superscript 1), 

where I = 2,3,4. This is comprised of three independent boson Hamiltonian corresponding 
to "FM" spin, orbital, and spin-orbital waves. 

For the AF SU(4) system, the classical ground state may be obtained by replac- 
ing the operators in (j3|) by their expectation values with respect to states of the form 

4> T ^j ; and then minimizing the energy. Since 
Xmm (Sn{i)) (Sm(j)) > 0, the classical minimum for each bond is obtained by having the 
two sites connected satisfy (<pi\ \<j>j) = 0. This can be achieved by having any two sites 
connected by a bond having different flavors. Assuming this can be done for all bonds 
(unfrustrated) , the classical ground state is then identical to that of the AF 4-state Pott's 
model. For concreteness, let us consider for now the square lattice with n.n. coupling J and 
next n.n. coupling J', with J > 2J' . Out of the degenerate manifold of classical ground 
states, we select for illustration the one obtained by dividing the lattice into 4 interpenetrat- 
ing square sublattices and having all sites on sublattice a be in the flavor state \oc) (see Fig. 
la), i.e.(S™(i)) = S nm 5 nai , where a* denotes the sublattice the site i is on. Expanding @ 
to the leading order in M, we have H = J2mn H mn , where mn are the pairs of the 4-states 
in (|1|) with m 7^ n, and 

H mn = £ Mb^ii^ii) + b^(j)bl(j) + KH*WSj) + h.c.]}, (5) 

hi 

where % and j are summed over all the sites in sublattices m and n respectively. Note that 
the different pairs of H mn are decoupled as a consequence of the SU (4) symmetry. A boson 
b™ at sublattice m only coupled to its "mirror" boson 6^ at sublattice n. This allows us 
to solve each H mn separately, and H mn is a simple Bogoliubov problem identical to the 
usual spin wave theory in the HAH. The spin (s), orbital (r), and spin-orbital (sr) wave 
dispersions are given by oo s (k) = 2J^Jl — cos 2 k y , uo T (k) = 2Jy/l — cos 2 k x , and oj ST (k) = 
4 J' 'Jl — cos 2 k x cos 2 k y . Note that the spin and the orbital excitations are ID in spite of the 
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I*) = Uil^i, so that (Si 



underlying translationally invariant 2D ordering pattern. The ID density of states can be 
understood from the ordering pattern of the classical ground state in Fig. la. There is no 
energy cost if we take any vertical or horizontal line and interchange the flavors on the two 
sublattices (creating a line defect). The disordering effect due to quantum fluctuations of 
these ID like excitations should act just like those in ID quantum spin chains and destroy 
the LRO. Of course, a linearlized theory can only be used as a guide to the true situation. 



l4[13H Nevertheless the result here is consistent with the previous results that the ground 
state of the AF SU(4) model is a flavor liquid on many 2D lattices. In the limit J' — > 0, 
point defects ( e.g., flip flavor from 1 to 4 on a site) can be created with no cost in energy, 
thus oj st — > 0. Note that for J' > 0, the classical ground state does not have finite entropy 
per site in the thermodynamic limit, but the classical order is destroyed by inifnitesimal 
quantum fluctuations (all M < oo) . 

While the details of the above depend on the particular ordered state we choose, the 
central results of no mixing between different waves and existence of the excitations with 
reduced dimensionality are consequences of the SU (4) symmetry. We now study the effects 
of deviation from this symmetry. We will illustrate the physics using an example which is of 
particular interest, the 2D honeycomb lattice, corresponding to a single plane of V 2 3 , whose 
magnetic behavior Castellani et. al. || sought to explain by invoking orbital degeneracy. 
The experimentally determined magnetic structure of the ordered state is consistent with 
that of 4-sublattice ordering discussed above with the lattice appropriately modified ||, as 
shown in Fig. lb. The full Hamiltonian when Hund's rule and anisotropy are taken into 
account is quite complicated. Since we are principally interested in how the physics changes 
away from the <S77(4) symmetry in ways that are independent of the precise model, we will 
ignore Hund's rule and consider a Hamiltonian including orbital anisotropy, 

where (ij) denotes n.n. pairs, m,n,m' ,n' are summed from 1 to 4, and is the hopping 



integral from state | m) at site i to state | n) at site j [pj| . Because of orbital anisotropy, 



t^] is diagonal in spin space, but not in orbital space. The hopping matrix for different 
orientational bonds (i'j') and (ij) are related || by a rotational transformation. Note that 
while ^ can be made diagonal on any one bond, it cannot be simulataneously done for all 
bonds. We choose the atomic orbitals with the hopping matrix diagonal for the "horizontal" 
bonds with eigenvalues t < and t > . The quantity r\ = t > ~*^ is a measure of the strength of 
the anisotropy. Expansion of the Hamiltonian in M leads to the flavor wave Hamiltonian, 

H % = jjY. E 1^(0 W)*S£,*i£} + E (b^(i)b a /\j)t^l + h.c.)\. (?) 

Consider first the S77(4) limit, = £<5 m , n , namely t] = 0. is reduced to the form 
of (H) with J = 2t 2 /U, and the spin, orbital, and spin-orbital waves are decoupled. The 
different connectivity of the honeycomb lattice and the ordering pattern in Fig. lb results 
in some modifications from that in the sqaure lattice in the dispersion of the excitations. 
Flavors 1 and 4 (or 2 and 3) are connected as zig-zagging vertical chains, and spin-orbital 
waves are ID like with dispersion cu ST (k) = J sin(v / 3|^ y |/2). Flavors 1 and 2 (or 3 and 4) are 
unconnected, and spin waves are localized on a site and have zero excitation energy, uj s = 0. 
Flavors 2 and 4 (or 1 and 3) are connected as two-site pairs which are decoupled from other 
pairs. The orbital modes (bf, bf) are thus also localized with uj t = 0. 

With anisotropy, the situation changes significantly. Spin is still conserved, but orbital 
no longer. Thus, spin modes (62, and spin-orbital modes (b\, b^) are now mixed, and 
therefore both neutron scattering active. However, for convenience, we continue to label 
them as spin and spin-orbital modes. Orbital information can thus be obtained indirectly 
from neutron scattering experiment. For H hc , the dispersion of the spin-orbital mode remains 
ID like. The lack of coupling between zig-zagging vertical chains can be understood simply 
as follows. Because the hopping matrix is diagonal on the horizontal bonds, the coupling 
between the two sites connected by such a bond in (0) contains only terms of the form 

(i) S™(j). Thus, this coupling conserves S%(i)+S%(j), and hence the chains are decoupled 
for spin-orbital and also spin modes for the same reason as the SU(4) case. For spin wave 
mode, the dispersion in fact remains dispersionless with zero excitation energy, indicating 
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they remain localized, albeit no longer on a single site. Pure orbital excitations remain 
localized on the two sites of each horizontal bond and decoupled from others. 

In Fig. 2a we show the energy dispersion as a function of k y for the spin and spin- 
orbital modes, which are the modes that can be probed by neutron scattering. Excitation 
energy is shown in units of J = 2t>/£7. There is no dispersion with k x since the chains are 
decoupled for these modes. The four purely spin modes {b\, b\, b\, bf) in the S , f/(4) limit 
remains dispersionless with anisotropy but with some spin-orbital characteristics mixed in. 
The four spin-orbital modes^, 6§, b%, bf), which for SU(A) were all degenerate, are split 
with anisotropy into two two-fold degenerate branches with the upper branch no longer a 
Goldstone mode as k y — > 0. If the effects of Hund's rule are included with effective strength 
Jh, chains will be coupled, and all the modes will have 2D dispersion. Furthermore, any 
accidental degeneracy and/or zero energy excitations will be removed. The spectrum will 
consist of two degenerate Goldstone spin waves whose dispersion is linear for small k with 
a slope whose scale is determined by Jh', two other degenerate spin waves with a gap and 
dispersion also detemined by Jh', two degenerate spin-orbital waves with a gap and dispersion 
in k x of order J H , but whose dispersion in k y is of order J; two other degenerate spin-orbital 
modes with gap of order Ji], and k x and k y dispersion of order J H and J respectively. 
The twofold degeneracy in each case is required by the remaining SU(2) spin symmetry 
of the Hamitonian. These features are qualitatively in agreement with inelastic neutron 
scattering data on V 2 3 which shows that the effective in-plane spin-spin coupling between 
parallel spins to be considerably weaker than that between antiparallel ones. Details of the 
calculation including Hund's rule and interplane coupling of V 2 3 will be discussed in a later 
publication. 

Because of the mixing between spin and spin-orbital waves, neutron scattering exper- 
iments will couple to all the modes mentioned above. The intensity of inelastic neutron 
scattering cross section is proportional to Jm^u), where x is the transeverse susceptibil- 
ity. Within the flavor wave approximation, the spin lowering operators S~(i) = 82(1)+ S%(i), 
which is if « £ sublattice 1 or 3, and if i G 2 or 4, and similarly for the raising 
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operator. The relative intensity can thus be calculated straightforwardly using the Bogoli- 
ubov transformation. We show an example of this for T = in Fig. 2b for i] = 0.43 and 
without Hund's rule. For finite T, the intensity of each branch needs to be multiplied by 

(l + 2/(e^-l)). 

In summary, the Holstein Primakoff transformation has been generalized to develop a 
quantum flavor wave theory for spin systems with orbital degeneracy. In addition to spin 
and orbital bosonic excitations, proximity to SU (4) symmetry gives raie to a third mode of 
boson excitations observable by neutron scattering. 

We thank W. Bao and I. Affelck for many interesting discussions. The work is in part 
supported by DOE grant DE-FG03-98ER45687. 
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FIGURES 

FIG. 1. Four-sublattice classical ground states considered in the text, a) Square lattice; b) 
Honeycomb lattice. 

FIG. 2. a) Spin (triangules) and spin-orbital (circles and squares) wave spectra with orbital 
anisotropy r\ = 0.43 for the honeycomb lattice as a function of k for wavevector (k x ,k). Dashed 
curve shows the (4-fold degenerate) spin-orbital excitations in the SU(4) limit. b) Calculated 
neutron scattering intensities (arbitray units) for the spin (triangles) and spin-orbital (circles and 
squares) modes 
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